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Abstract 

We image spatial distributions of Xe 9+ ions in the focus of a laser beam of ultrashort, intense 
pulses in all three dimensions, with a resolution of ~3 /jm and ~12 urn in the two transverse 
directions. This allows for studying ionization processes without spatially averaging ion yields. 
Our in situ ion imaging is also useful to analyze focal intensity profiles and to investigate the 
transverse modal purity of tightly focused beams of complex light. As an example, the intensity 
profile of a Hermite- Gaussian beam mode HG^o recorded with ions is found to be in good agreement 
with optical images. 
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After the development of chirped pulse amplification [l| in the late 1980s and the discov- 
ery of self-mode-locking[2] in 1991, laser sources emitting femtosecond-duration pulses have 
become routinely available. Laboratories around the world have used such sources for exten- 
sive studies of intense-field processes such as multiphoton ionization, tunneling ionization, 
and above-threshold ionization of atoms and of molecules |3j]. The required atomic electric 
field of about 10 9 V/cm or more is created by focusing the ultrashort pulses in a target 
gas. Currently, much attention is geared toward the experimental realization and character- 
ization of attosecond pulses using intense-field processes 0, 0]. With such extremely short 
pulses the dynamics of electrons in atoms and molecules can be followed on their own natural 

n 

timescale; this now emerging field is called 'attophysics'jgj. 

A direct consequence of carrying out intense-field experiments with focused laser light is 
the presence of a broad range of peak intensities across the focal region. To understand some 
of the experimental consequences of this let us consider a process for which the probability 
to produce a certain product particle is P(I) when the pulse has peak intensity I. Following 
common practice, we might study the process by adjusting the laser pulse energy such that 
the highest peak intensity found in the focus equals Jo and then collecting product particles 
from all over the focal region. Doing so, we obtain a volume-integrated product yield Y(I ) 
given by 
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F(I )oc !° P{I) 
Jo 



in which V(I,Iq) is the volume in which the local peak intensity exceeds some value /Q]. 
Because of the volumetric weighting in Eq. (pQ), the measured yield Y(Iq) differs from the 
probability P(I) under investigation. For instance, when I is increased, Y(I ) ultimately be- 
comes proportional to I^ 2 once the saturation intensity (for which P(Iq) — > 1) is exceeded 8]. 
This masks the decrease of P(I) that is often expected for intensities greater than the sat- 
uration intensity due to competing higher-order processes. Also, subtle features in P(I) 
0, [h3| tend to be washed out by volume integration. 

In this Letter, we present spatially resolved images of 3-D densities of ion charge states 
as they are created in focused ultrashort pulses. Thus, what we have realized for use in 
intense-field ionization research is a photodynamical test tube, i.e. a volume of known, jum- 
sized dimensions exposed to peak intensities that are essentially constant across the volume. 
This allows measuring P(I) without spatial integration: one simply defines a sufficiently 
limited spatial region of interest in the ion image and determines the yields originating from 



this region. 

Another goal of our investigations is to analyze focal intensity distributions. These can 
be inferred because volume elements having the same peak intensity contain identical ion 
charge state distributions. The target gas then acts as an in situ sensor of local intensity. 
In this Letter we analyze images, made with Xe ions, of the transverse profile of an intense, 
ultrashort focused Hermite- Gaussian HGo,i transverse beam mode. The HG modes are 
related to the Laguerre-Gaussian modes LG Pj £ll|. We are particularly interested in LG 



modes 
(OAM) 



because they carry a sharply defined amount of optical orbital angular momentum 



121 ]. and we wish to study the effect of this quantity on intense- field processes. For 



this, we need modally pure LG modes. In situ ion imaging will be helpful to assess the 
degree of modal purity: modal interferences due to contaminating mode components (which 
have different Gouy phases [ill]) should show up in the ion image. 

Other techniques have been reported to image foci and/or circumvent volumetric weight- 
ing. First, the intensity-selective scanning (ISS) technique @,Q,Q,QQ utilizes a narrow 
slit perpendicular to the laser propagation direction, so ion yields are still integrated over 
two dimensions. To obtain yields free of volume integration raw ISS data must be decon- 
volved, which adds noise and only works for Gaussian intensity distributions. Second, the 

, [isl has been shown to work only for a two- 
dimensional Gaussian profile, and is inherently prone to statistical noise. Lastly, a number 



of time-of-fli 



times 




it techniques allow ions created at different locations to arrive at different 



21]. The technique of [19[] was demonstrated to work well for ~200 /zm 
resolution [20|, but more systematic investigations are required[jjj to explore its applicability 
for resolutions of a few micrometer. The setup of 21[ only achieved a spatial resolution of 
0.68 mm, and requires cold targets. A focus with a waist of a few mm was dissected into 12 
adjacent detection volumes of a few 100 /zm. 9| 

Compared to the previous methods, the present method achieves a better spatial resolu- 
tion (~3 fim and ~12 /zm in the two transverse directions of the beam). It does not require 
a priori knowledge about the intensity distribution, and is not limited to Gaussian profiles. 
We show the first spatially resolved ion distributions in three dimensions with resolution of 
a few micrometer. 

A schematic of our setup is shown in Fig. [TJ We used an 800-nm Ti:sapphire laser system 
(Spectra Physics Tsunami/Spitfire) delivering ~50-fs pulses at a 1-kHz repetition rate with 
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FIG. 1: Schematic of our setup (see text for details). 

an average output power of about 1.2 W. We focus the radiation into an evacuated ionization 
chamber using an AR-coated UV-grade fused silica plano-convex lens with a nominal focal 
length of 22.7 cm for 800 nm. We lowered the background pressure in the ionization chamber 
to below 10~ 9 mbar using a liquid-N 2 trap. We then admitted Xe gas (purity 99.999%) to 
a total pressure of ~ 5.0 x 10~ 7 mbar. To detect the Xe ions created in the focus, we use 
a reflectron-type tof ion mass spectrometer. The focus is positioned between two parallel 
plates. One of these is at ground potential and contains the spectrometer entrance slit; the 
other plate ('repeller') carries a positive voltage Va of about +1.5 kV and is at a distance 
d of about 3 mm away. The uniform electrostatic field between the plates causes target gas 
cations created in the focus to be accelerated toward and pass through the slit. The slit 
is rectangular, measuring (see coordinate system in Fig. [1]) Az ~400 /iin by Ay ~12 fim. 
Assuming their initial kinetic energy is negligible, ions created at a distance x away from 
the slit enter the tof tube with kinetic energy proportional to (x/cI)Vr. The total tof of 
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an ion with mass m = Mu (u is 1 a.m.u. ~ 1.67 x 10 27 kg) and charge q = Qe (e is the 
elementary charge ~ 1.60 x 10~ 19 C) equals 




i = ^r W . (2) 

In this equation, r(x) is a reference tof of an ion with M = Q (e.g. H + ) that was created 
at position x; this time depends on the (adjustable) distance d and the fixed distances d\, 
d 2 , and £ = £\ + £2 between the meshes M 1; M 2 , M 3 and the slit, and on the electrostatic 
potentials V2 and V3 we apply to the meshes (see Fig. [Q note that V\ = V). The ion 
trajectories inside the reflectron have a turning point, and the ions are detected on a multi- 
channel-plate ion detector (mcp in Fig. [I]). To improve the detection efficiency the ions are 
post-accelerated by applying a potential Vp OSt ~ 8 kV over a distance d post before they hit 
the mcp. The mcp was a Galileo chevron type with a time resolution of At f» 1.5 ns. To 
record the ion signal from the mcp we used a FAST ComTec model P7886 2-GHz counting 
card with time bins of 500 ps. 

To reconstruct three-dimensional ion distributions, we note that the entrance slit already 
clips the acceptance volume in the two dimensions y and z to the sizes Ay and Az of 
the slit. We reconstruct the initial positions of ions in the third dimension (x) from their 
tof spectrum (Eq. [2]). To obtain a one-to-one correspondence between initial location x 
of an ion and its tof for each M/Q ratio, we adjust the mesh potentials V2 and V3 so 
that a monotonically increasing r(x) results. For a fixed accumulation time, the signal Sj 
recorded in time bin j is proportional to the number of ions with tof between tj = jAt and 
tj + i = (j + 1) At (At is time bin width). These ions must all originate from a small spatial 
interval Axj with tj = yjM/Q r{xj) and t j+t = y/M/Q r(x j+1 ). The density 

of ions Nj created in the spatial interval Axj is thus 

Nj oc =^Sj oc (3) 
J Axj 3 Axj K ' 

The second proportionality follows because all Atj are equal. We may think of the recorded 
tof spectrum Sj as the image in time of the density Nj of ions in space. The spatial 
resolution (in the x direction) is 

Ax-^M- ^ (4) 
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FIG. 2: (Color online) Measured distribution of 132 Xe ions in an xy plane through the focus, with 
Xe + in red, Xe 2+ in green, Xe 3+ in blue, and Xe 4+ in red, as indicated. The isotope 132 was used 
because it is the most abundant. The signal on the left is due to 131 Xe + . Cross sections along the 
dashed line are shown on the right. The observed nesting of the charge states occurs because the 
peak intensity is higher at locations closer to the center of the focus. A charge state that dominates 
at lower-intensity locations is depleted in favor of the next-higher charge state at higher-intensity 
locations. 

For the typical settings we used, the derivative dr/dx smoothly drops from about 0.250 
ns//xm (at x = 1.9 mm) to about 0.200 ns//xm (at x = 2.0 mm). To measure r(x) we 
used a focused attenuated beam, which generates only Xe + . We positioned this focus at 
various x positions using a translation stage with /im precision, and recorded the tof of 
the Xe + . The obtained experimental t(x) matched the theoretically predicted curve. Using 
Eq. (j3J) with a time resolution At = 1.5 ns as set by our mcp, we find our theoretical spatial 
resolution for 132 Xe + ranges between 0.52 /im (at 1.9 mm) and 0.66 /im (at 2.0 mm). For 
132 Xe 7+ , these values would be 1.4 /im and 1.7 /im, respectively. A more realistic estimation 
of the spatial resolution accounts for the thermal motion of the Xe gas, the dimensions of 
the spectrometer entrance slit, mechanical vibrations, and the frequency bandwidth of our 
detection electronics. A first analysis of these factors indicates the spatial resolution in the 
x-direction is about 3 /im. In the y and z directions the resolutions are mainly determined 
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FIG. 3: (Color online) Measured distribution of 132 Xe ions in the xz plane, with Xe + in red, Xe 2+ 
in green, Xe 3+ in blue, and Xe 4+ in red, as indicated. The dashed line is the optical axis (beam 
propagates to the right). Note the difference in scale between the x and z axes. The small angle of 
23 arcmin between the optical axis and the z axis is due to a slight misalignment of our setup. The 
signal at the bottom of the image is due to 131 Xe ions. The observed nesting of the charge states 
occurs because the peak intensity is higher at locations closer to the center of the focus. A charge 
state that dominates at lower-intensity locations is depleted in favor of the next-higher charge state 
at higher-intensity locations. To avoid overlap between adjacent isotopes we used slightly different 
voltages than those used to obtain Fig. [2] 

by the slit dimensions, so 12 fim and 400 /im, respectively. Slits narrower than 12 /im could 
be used, but the resolution in the y dimension is ultimately limited by thermal effects, in 
particular for lower charge states. We chose the 400-/im slit length as a compromise between 
count rate and the need to keep the slit shorter than the Rayleigh range of our focused beam, 
which is a few millimeter. Pressure-dependent measurements indicated that space charge 
effects are absent for the pressure we used (~ 5.0 x 10 -7 mbar). 
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Figures [2] and [3] show 132 Xe ion distributions in the focus reconstructed using Eqs. [3] and 
HI The different charge states are color-coded: Xe + (red), Xe 2+ (green), Xe 3+ (blue), and 
Xe 4+ (red). The image in Fig. [2] is a cross section in the xy plane through the focal region, 
obtained by moving the focus up (y direction) in steps of 10 /im using a stepper motor, and 
recording a tof spectrum for each height. The image in Fig. [3] is a cross section in the 
xz plane, obtained by moving the focus along the z direction in steps of 200 fim. Setting 
up the experiment, we noticed that visible imperfections in the unfocused beam profile go 
hand in hand with irregular features in the ion charge state distributions. Using this as 
a guideline, we improved the mode quality of the laser by reducing the aperture of an iris 
in the regenerative cavity of the amplifier, at the expense of ~0.1 W. The nesting of the 
various charge states in Figs. [2] and [3] indicates the competition between the processes that 
create them. The four line graphs in Fig. [2] show the ion yields measured along the dashed 
line in Fig. [2] (transverse to the laser propagation, same scale). If one moves toward the 
center of the focus, the yield of each charge state Xe q+ rises, saturates, and is then depleted 
in favor of the next charge state Xe < - <?+1 - )+ in a sequential way, so that a nested pattern of 
charge states results. The depletion shown in Fig. [2] is in sharp contrast with the 3/2-power 
rise observed in traditional methods. From the data shown in the line graphs, we measured 
that the count rate of Xe + and Xe 2+ in the center of the focus is only a few percent of 
that of the peaks. Such a degree of observed depletion has no precedence in literature, and 
underlines our ability to record P(I) free of volumetric weighting. 

To demonstrate that our ion imaging technique is not limited to only Gaussian intensity 
profiles, we also imaged a focused beam of ultrashort pulses in the Hermite-Gaussian HGi n 
transverse mode[rj| (Fig. H]). We created this mode using holographic techniques [221. |23 |. 
The intensity of the HG^o mode is 

I(x, y, z) oc — — ^ exp — - , 5 

in which w(z) 2 = Wq (1 + z 2 /zq) describes the self-similar transverse scaling of the intensity 
profile (w is the beam waist and z is the Rayleigh range) . The profile consists of two 
lobes (see inset a in Fig. H]). For y = (dashed line in inset) the maxima of these lobes are 
separated in the x direction by a distance 



Ax{z) = Ax Jl + z 2 /z 2 (6) 
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(see inset 6), with Ax the separation at the waist (z=0). Again using Eqs. [3] and IU we 
reconstructed Ax(z) (circles in Fig. H]) from tof spectra. 

To verify that this reconstruction is correct we also focused the beam in air using the 
same lens as in the ionization experiment, and recorded the beam profile as a function 
of z optically (triangles in Fig. H]) using a CCD camera. In contrast to the ion imaging 
experiments — taken at full power — this optical experiment required substantial attenuation, 
to protect the camera from damage and to avoid nonlinear optical effects in air. The ion 
image and optical image agree within the error limits, which we estimate to be ~8.4 fim 
(pixel size) for the CCD data and ~3 /im for the ion data. This agreement allows a fit to 
Eq. flS]) (solid curve in Fig. 0J. 

In summary, we used a time-of-flight technique to image ion distributions in the focus of 
a laser beam of ultrafast pulses. We showed that this (a) allows measurements of ion yields 
without volumetric weighting, and (b) facilitates investigations of focal intensity distribu- 
tions, with the target atoms acting as in situ intensity sensors. 

This material is based upon work supported by the National Science Foundation under 
Grant No. PHY-0355235. The authors gratefully acknowledge the Max-Planck-Institute of 
Quantum Optics in Germany, in particular Dr. Hartmut Schroder, for generously lending us 
the reflectron used in this work. 
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